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Introduction 

In relativistic kinetic theory, global dynamics of several kinds of charged 
and uncharged matter remain an active research domain in General Relativity 
(GR), in which cosmology plays one of the central roles, by coupling various 
matter fields to the Einstein equations to provide mathematical explanations 
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in response to new astrophysical observations. In this context, spatially ho- 
mogeneous phenomena such as the one we consider in the present paper are 
relevant. There are several reasons why it is of interest to consider the Einstein 
equations with cosmological constant and to couple the equations to a massive 
scalar field. 

• Astrophysical observations have made evident the fact that, even in the 
presence of material bodies, the gravitational field can propagate through space 
at the speed of the light, analogously to electromagnetic waves. A mathema- 
tical way to model this phenomenon is to couple a scalar field to the Einstein 
equations. Let us recall that the Nobel prize of Physics 1993 was awarded for 
works on this subject. More details on this question can be found in [5], [21]. 
In fact, several authors realized the interest of coupling scalar field to other 
fields equations; see for instance [i], [18], [20], [24], [8]. 

• Now our motivation for considering the Einstein equations with a cosmo- 
logical constant A is due to the fact that astrophysical observations, based on 
luminosity via redshift plots of some far away objets such as Supernovae-Ia, 
have made evident the fact that the expansion of the universe is accelerating, 
as foreseen by E.P. Hubble. A classical mathematical tool to model this phe- 
nomenon is to include the cosmological constant A in the Einstein Equations. 
Several authors did it in the case A > 0; see for instance [14], [7], [15], [22], 
1 23], [26]. In the present paper, we prove that, in the presence of a massive 
scalar field, this result can be extended, not only to the case A = 0, but also to 
the case A > —a 2 , where a > is a constant depending only on the potential 
of the massive scalar field. This result extends and completes those of [20]. 

Also recall that the recent Nobel prize of Physics, 2011, was awarded to 
three Astrophysicists for their advanced research on this phenomenon of ac- 
celerated expansion of the universe. 

In fact, we must point out that, the notion of "dark energy" was intro- 
duced in order to provide a physical explanation to this phenomenon, but the 
physical structure of this hypothetical form of energy which is unknown in the 
laboratories remains an open question in modern cosmology; so is the question 
of "dark matter". Also notice that the scalar fields are considered to be a 
mechanism producing accelerated models, non only in "inflation", which is a 
variant of the Big-Bang theory including now a very short period of very high 
acceleration, but also in the primordial universe. 

In this paper, we consider the 3+1 formulation of the Einstein equations, 
which allows to interpret the fields equations as the time history of the first 
and second fundamental forms of the 3-hypersurfaces of constant times slices, 
foliating the space-time. The background space-time is any Bianchi space-time 
type I to VIII, since it is proved in [ ] that Bianchi IX such as the Kantowski- 
Sachs space-time, develops curvature singularities in a finite proper time, and 
this constitutes a major obstacle to our goal which is to proved global in time 
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existence of solutions. We prove in this paper that if the initial value of the 
mean curvature is strictly negative and if A > —a 2 , then the coupled Einstein- 
Maxwell-Scalar Field System has a global in time-solution. We investigate the 
asymptotic behaviour which reveals an exponential growth of the gravitational 
potentials, confirming the accelerated expansion of the universe. We prove the 
geodesic completeness and we were able to show that the considered model 
always satisfies the weak and the dominant energy conditions; we prove that, 
if A > f3 2 where > is a constant depending only on the mean curvature 
of the space-time, then the considered model also satisfies the strong energy 
condition. 

The paper is organized as follows: 

• In section 1, we introduce the coupled system and we give some preli- 
minary results. 

• In section 2, we study the constraints equations, the mean curvature and 
we introduce the Cauchy problem. 

• In section 3, we prove the local and the global existence of solutions. 

• In section 4, we study the asymptotic behaviour. 

• In section 5, we study the geodesic completeness. 

• In section 6, we study the energy conditions. 

• Section 7 is the appendix to which we refer for the details of the proofs 
of some important results. 

1 Equations and preliminary results 

• Unless otherwise specified, Greek indices a, (3, A, range from to 3 
and Latin indices i, j, k, from 1 to 3. We adopt the Einstein summa- 
tion convention a a b a = J2 a ab a . 

a 

We consider a time-oriented space-time (M, g) where M is a four- 
dimensional manifold and g the metric tensor of lorentzian signature 
(—,+,+,+). The model adopted for our study is any Bianchi space- 
time from type I to VIII, and, following [17], [19], [25], [27], who studied 
the question, we take M on the form: M = K x G, where G is a three 
dimensional simply connected Lie group. We take g on the form: 

g = -dt 2 + gij{ty <g> e j (1) 

where {ej} is a left invariant frame on G and {e 1 } is the dual frame; (p^) 
is a positive definite Riemannian 3-metric depending only on t. Now the 
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vector n = d t being orthogonal to G, we complete the frame (ei) on G, 
to obtain a frame (n, e 7 ) on M. We have 



g(n,ei) 

confirming the form (1) of g. 



goi 



0. 



The Einstein-Maxwell-Massive scalar field system with cosmological con- 
stant A, which rules the evolution of the considered charged pure matter 
can be written, following [ ], and denoting with a tilda (~) quantities on 
M: 

R a p - ]-Rg a p + Ag a p = 87r(T Q( g + r a/3 + pu a u p ) (2) 



V„F Q/3 = An J 13 



-yF, 



a/3 



(3) 
(4) 



where: 



(2) are the Einstein equations, basic equations in GR for the un- 
known metric tensor g = (<jf Qj g); R a p is the Ricci tensor, contracted 
of the curvature tensor, R = g a/3 R a p is the scalar curvature; T a p 
and r a p whose expressions are given below, are respectively, the 
tensor associated to a massive scalar field $ which is an unknown 
function of the time t, and r a p the Maxwell tensor associated to the 
electromagnetic field F. 

(3) and (4) are the two sets of Maxwell equations, basic equations of 
Electromagnetism, written in covariant form for the electromagnetic 
field F = {F° l , Fij) which is a closed unknown antisymmetric 2- 
form, depending only on the time t. F 0t and F^ are respectively its 
electric and magnetic parts. T a/3 and f a p are defined by: 



T, 



a/3 



T, 



a/3 



V^V^ - -g aP 
\g a pF^F 



\fi + F a \Fp A 



(5) 
(6) 



where in (5) as in (3) and (4), V stands for the covariant derivative, 
or the Levi-Civita connection in g, m > is a given constant called 
the mass of the scalar field notice that |m 2 $ 2 represents the 
potential associated to the scalar field $. 

In (2), pu a Uj3 is the tensor associated to the considered charged 
pure matter, with p > an unknown scalar function of the time t, 
standing for the pure matter proper density, and u = (ffl) a time- 
like future pointing unit vector which is an unknown function of the 
time t, representing the material velocities. 
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- In (3), J = (J 73 ) stands for the Maxwell current generated by the 
charged particles of pure matter and defined by: 

J? = eu p (7) 

in which e > is an unknown scalar function of the time t, standing 
for the proper charged density of the charged particles. 
Notice that the Maxwell equations (4) are just the covariant nota- 
tion of the relation dF = 0, since F is a closed 2-form. 
Now it is well known, see [12], that the electromagnetic field de- 
viates the trajectories of the charged particles which are no longer 
the geodesies of the space-time as in the empty case, but the solu- 
tions of the following differential system of current flow: 

u a V a u p = 4tt^F a V (8) 

which reduces to the usual geodesies system when F — 0. 

• To study the Einstein equations (2), we adopt the 3+1 formulation, 
which allows to interpret these equations as the evolution in time of the 
triplet (T,t,9t,k t ), where £t = {t} x G, g t = (fi^(t)) stands for the 
first fundamental form induced on J2t by g, k t = (&y(t)) is the second 
fundamental form defined in the present case by: 

hj = --dtga (9) 

where we adopt the ADM/MTV convention see [1], [6]. 
By the 3+1 formulation, the Einstein equations (2) which originally, in 
the homogeneous case we consider, are a non linear second order differen- 
tial system in gij, can be written as an equivalent first order differential 
system in (g^, ky) to which standard theory applies. 
We now introduce a quantity which will play a central role, namely, the 
mean curvature of the space-time which is a scalar function denoted by 
H and defined by: 

H = g ij k i3 (10) 

that is the trace of (fcy). Next, since u = (u^) is a unit vector which 
means g(u, u) = u a u a = — 1, we deduce from the expression (1) of g that: 

u° = y/l + giftv? (11) 

because u is future pointing, and (11) shows that g^ and u % determine 
u°, and, since (g^) is positive definite, that u° > 1. 
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It is important to note that (e,) is not always the natural frame (<9j) on 
G. We set (e ) = (n, e*) = (dt,ei): that is eo = <9 t and = e^. The 
frames (e a ) and (<9 a ) are then link by: 



e X a dx (12) 



where 



S° = l; e« = ?o = 0; e*=ej. (13) 

Recall that the structure constants of the Lie algebra q of the Lie 
group G are defined by: 

fe, e,] = Cje* (14) 

where [ , ] denotes the Lie brackets of g. Since the Lie brackets are 
antisymmetric, C\a is antisymmetric with respect to i, j that is: 

C% = -C*. (15) 

Now the Ricci rotation coefficients 7^ associated to the Levi-Civita con- 
nection V are defined by: 

= lapex- (16) 

We set 

[e ai ep] = C^e x (17) 

where: 

U aB — U 0B — u a0 — U ' U ij ~ 

The general expression of 7^ in term of C^b an d 9aB can be find in [ 
or in [3] p. 301. It will be enough to extract and mention here only those 
of the Ricci rotation coefficients used in the present paper, namely: 

7oo = 7°o = 7oi = °; % = % = iL = ~ki 



lij ■ lij 2^ 



(19) 



Notice that the 7^ are not to be confused with the usual Christoffel 
symbols r** associated to the natural frame (d\). For instance, (19) 
gives, using Cg = -Cg: 

7i = c&; i\i-i l n = G \i (20) 

which shows that, at the contrary of the T^, the 7^ are not symmetric 
with respect to a and (3. We have the following formulae, analogous to 
the natural frame case: 




e a {T B ) - 7*,T A . 



a x (21) 
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Also mention the useful formula we prove in Appendix A\\ 

d 
dt 

where V is the Levi-Civita connection on (G,g). 



(llj) = V'fcy - - V<*j (22) 



A direct calculation using (21) shows that the components of the tensor 
T aj 3 defined by (5) are given, setting $ = by: 

Too = \& + l -m 2 <$> 2 - f 0l = 0; % = l -g^ 2 - m 2 $ 2 ). (23) 

Concerning the Maxwell tensor r a p defined by (6), we first obtain by a 
direct calculation: 



f F X »F^ = -2 gij F^FV + g ik g^F iL F M ; 
Fi\F, x = -g ikgjl F ok F 01 + g kl F ik F jf , 



(24) 



Fn\F- x = -F uF ok - 
F QX F x = g l3 F *F°i; 

from where we deduce, using (6), that: 

foo = \gaF<*F® + ^g ik g jl F k iFij 
Toj = —F ok Fjk 

Tij = {\gii9u ~ 9 i kg J i)F 0k F 01 - \g l0 g km g nl F kl F mn + g kl F tk F 3l . (25) 

Next, to derive the equations for the scalar field $ and the matter density 
p, we use the conservation laws: 

V Q f aP + V Q f a/3 + V a (pu a uP) = 0. (26) 

A direct calculation using (5), (6), (21) and the Maxwell equation (4) 
gives: 

V Q f Q/3 = V^$(V a V Q $ -m 2 $); V a f aP = F p x V a F aX (27) 

where V a V a often denoted □- is the d'Alembertian or the wave operator. 
Now (26) and (27) give, using the Maxwell equation (3): 

V^$(V a V a $ - m 2 $) + 4<KeF x u x + u p V a (pu a ) + (pu a )V a u^ = (28) 

(28), reduces using (8) to: 

V^(V Q V a $ - m 2 $) + u fi V a (pu a ) = 0. (29) 
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But it is easily seen that V*$ = 0; so (29) gives for (3 = i, u l V a (pu a ) = 
which is satisfied for every i — 1, 2, 3 if: 

w a ( P u a ) = o. (30) 

Now (29) gives, since V J $ = and using (30): 

V°$(V Q V a $ - m 2 $) = 0. (31) 

It then appears that the conservation laws (26) will be satisfied if equa- 
tions (30) and (31) in p and <3> are. 

Now a direct calculation shows that equation (30) in p writes: 



p(0)5°(0) - /•' »' : 
P = =n exp 



which solves at once over [0,t], t > 0, to give: 

^i-qA)(*H- ( 32 ) 

(32) shows that (p(0) > 0) =>■ (p > 0). In what follows we set: 

p(0) > 0. (33) 
Next, it is easily seen, using (21), that 

V a V Q $ = -$ + H$, 
so that equation (31) in $ can be written, using V°$ = ^ 00 Vo$ = — $: 

$$ _ #($) 2 + m 2 $$ = 0. (34) 
To study this non linear second order equation in $, we set: 

U=\m 2 - (35) 

We choose to look for a non-decreasing scalar field which means $ > 0; 
(35) then gives: 

6 = V2U. (36) 

Next, to derive equation in e, we use the Maxwell equation (3), which 
gives, using the identity V^V^F ^ = and the expression (7) of 

V a (eu a ) = (37) 



Einstein-Maxwell-Massive Scalar Field System in 3+1 formulation 



(37), yields, similarly to (30): 

'•' ■ v? 



e(0)5°(0) 

exp 



u° 



ft n,3 

\ (K ~ ClX)(s)ds (38 



by (38): (e(0) > 0) e > 0. In what follows we set: 

e(0) > 0. (39) 

We now write in details the Maxwell equations (3), (4) in F; (3) gives 
for p = i and (3 = 0, and using the expression (7) of J 13 

f V a F ai = Aneu 1 (40) 
\ V a F a0 = Aireu . (41) 

- Applying formulae (21) yields: 

V a F m = i a F ai + XxF M + laxF aX - (42) 

But by a direct calculation using (19), (20) the properties 
C* = -Cg, F ij = -F ji and since 

e a F ai = d F 0i = F 0i 
we deduce from (40) and (42), the equation in F 0t : 

poi = H foi _ C j k pki _ ho) k F jk + Aiterf . (43) 

Now concerning (41), we have, using (19) and (20): 

V a F a0 = Cl k F k0 , so that (41) gives the constraints equation: 

Ci k F ok + Aneu = 0. (44) 

- Next, observe that the Maxwell equations (4) split into the two sets: 

VoFij + ViFjQ + VjF 0l = 0; ViF jk + VjF ki + VkFj = (45) 
the first equation (45) gives the equation in F^: 

% = -C k F ok (46) 
and the second equation (45) gives the constraint equation: 

C\ 3 F kl + C l kl F 3l + C] k F a := C l {ij F k]l = (47) 
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Since all the indices are now fixed, we set from now, to simplify notations: 



?0i 



Fiji T 00 



TOO — TOO] T 0i — T 0i'i T ij 



Too', T 0i 





Tm ; XL- 



'-Oi 



T- 



(4? 



Finally, since the system in (u l ) is given by (8) for /? = i, it remains to 
explicit the Einstein equations in fcjj); the uncharged case F = is 
given by classical equations. The equations in the charged case F / 
were set up by [16], and we easily adapt to the present case. We are then 
led to the following evolution system in (g^, k^, E\ Fij,u l , $, U, p, e): 

(49) 



9ij 



2ki 



kij = Rij + Hkij - 2k]ku ~ 87r(Ty + r y + pu(Uj 
- Too - pu\ + g lm (Ti m + puiu m ) - Ag 



+47T 

E l = 

Fij — C^g k iE l 



'.i 



- [ q k g kl g im + -Cy l g km ) F lm + Aneu 1 



u 
U 



2ky 



i u j u k 



e „„• . e g d Fjiu j 



4-k-E 1 + 47T- 

p 



2HU-m z $V2U 

P = - 



p + 4ireg, 



k ■ 



MM 



j 



e + 47rg, 



E^v? e 2 



13 <v?y P 



(50) 

(51) 
(52) 

(53) 

(54) 
(55) 

(56) 
(57) 



where: 

- In (50) is the Ricci tensor associated to gij and given see [2] by: 

= 7z m 7j; - 7j7 7 mi - C mi 7 K ; (58) 

- In (53) Yjk 1S given by (19); 

- The equation (55) in U is given by (34) using the change of variable 
(35) which provides in the same time equation (54) in $ following 
the choice $ > 0; 

- Equation (53) in u % is given by (8) for (3 = i; 

- Equations (56) and (57) in p and e are given respectively by (30) 
and (37) in which w° is provided by (8) for (3 = 0, which gives: 
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Finally, (49) is provided by (9) and (50) is a direct adaptation of 
(14) in [16] to the present case. 



Next we have the following set of constraints equations: 

' R - kijk ij + H 2 = 167r(T 00 + roo + pul) + 2A 
V'kij = -87r(T j + r 0j + pu Uj) 
C[ijF k ]i := CljFu, + C l ki Fji + C l jk F u = 
{ C\ k E k + 4vreM° = 



(59) 
(60) 
(61) 
(62) 



where (61) and (62) are given by (47) and (44) whereas (59) (called the 
Hamiltonian constraint) in which R = g^Rij and (60) are easily deduced 
from (1.20) and (1.21) in [16]. 



2 Study of constraints and mean curvature: 
the Cauchy problem 

We first set up the evolution of the different quantities involved. We prove: 
Lemma 2.1. 

If the evolution system is satisfied, then we have: 



dH 

~dT ~ 
dH 2 

dt 



R + H 2 + kng i3 {T l3 + pUiUj) - 127r(T 00 + pu 2 ) - 87rr 00 - 3A. (63) 



2H 



R + H 2 -3A + ^g ij {T id + pu iUj ) - 127r(T 00 + pu 2 ) - 87rr 00 



j t {k ij kij) = 2k ij R, tj + 2H(k ij k ij - A) - I6n k ij (Tij + Tij + pu lUj ) 



(64) 
(65) 



d 
dt 

d 
dt 

d 
dt 



+87tH[-T 00 - pul + g lm (T lm + puiu m )}. 
(T 00 + Too + pul) = H(T 00 + r 00 + pu 2 ) + ^(Ty + r {j + pUiUj) 

-ViiT 01 + r 01 + pu°u l ). (66) 
(T 0j + r 0j + pu°u j ) = H(T 0j + r 0j + pu oUj ) + 2k{(T 0i + r 0i + pu'u ) 

-V,(T ij ' + r ij + pu i u j ). (67) 
R = 2k ij Rij - 2ViVik il . (68) 



Proof. 

See Appendix A 2 . 



□ 
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Lemma 2.2. 

Set- 



then 



Proof. 



' A = R — kijk* + H 2 - 167r(T 00 + too + pu 2 ) - 2A; 
Aj = V l kij + 87r(T 0j + T 0j + pu Uj); 
Aijk = C\,F k i + CtFji + CLF U ; 



B = C\ k E k + 4vreM C 



dA 

~dt 
dA l _ 

dt 

dAijk 

dt 
dB 

{ ~dt ~ 



2HA + 2gi m 1 k mj A k 



HA, 



= 
HB. 



(69) 
(70) 
(71) 
(72) 



1. ) (69) is a consequence of (64), (65), (66) and (68); 

2. ) (71) is a consequence of the evolution equation (52) in F^ and the Jacobi 

polynomial relation Ch^Cu = 0; 

3. ) To obtain (72), first set: 

= V a F af} - Averf (73) 

for (3 = i, given the evolution equation (40) in F° l = E l which is also 
written in the form (43) or (51), we have 



C l = 0. 



(74) 



Next, due to the identity V a VpF af) = and (37) i.e V^eu 13 ) = 0, we 
have 

= 0; (75) 
then by a direct calculation using (19), (21) and (74), (75) gives 

d t C° - HC° = (76) 

But (73) gives: 

C° = V a F a0 - Aneu = -C\ k F Qk - Aneu = -C\ k E k - Aixeu = -B. 
So (72) follows from (76). 
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4.) For the proof of (70) see Appendix A 3 . 

□ 

Proposition 2.3. 

The constraints (59), (60), (61) and (62) are satisfied in the whole domain of 
existence of the solutions of the evolution system (49), (50), (51), (52), (53), 
(54), (55), (56), (57), if and only if, they are satisfied for t = 0. 

Proof. 

Integration of (69), (70), (72) over [0,t], t > gives: 

Aj = Aj{0) exp( /' Hds); A tjk = A ijk (0); B = B(0) exp( / Hds)] 
Jo Jo 

hence (A^O) = A ljk (0) = B(0) = 0) <=^> (Aj = A ijk = B = 0). 
Now setting A k = in (69) gives A = 2HA which integrates over [0,t], t > 
to give A = A(0)exp(f*2Hds). Hence (A(0) = 0) (A = 0) and proposi- 
tion 2.3 follows. □ 

We suppose from now on that the constraints (59), (60), (61) and (62) are 
satisfied for t = 0. As consequence, we can use the constraints which can now 
be considered as properties of the solution of the evolution system. 

We now prove an important theorem on the mean curvature of the solutions 
of the evolution system. 

Theorem 2.4. 

Let $(0) > 0, A > -47rm 2 ($(0)) 2 be given, and suppose H(0) = (g ij kij)(0) < 
then H is uniformly bounded and we have: 

H(0) <H< - V / 3A + 127rm 2 ($(0)) 2 . (77) 

Proof. 

Denote by the traceless tensor associated to kij, i.e 



A direct calculation gives: 



H , . 

hj = -j9ij + cry- (78) 



fc% = l -H 2 + o-ija*. (79) 



We now use the evolution equation (63) in H, in which we use the Hamiltonian 
constraint (59) to express the quantity R + H 2 , (23) and (35) to express T 00 , 
Tij in terms of $ and U, and finally (79) to express kyk 1 ^ , to obtain: 

dH H 2 

— = — - A - 47rm 2 $ 2 + ayo*' + 16-kU + A^g^muj + ^pu\ + 87rr 00 . (80) 

QX o 



14 



JV. Noutchegueme and A. Nangue 



But since a^a^ > 0, g tj UiUj > 0, r 00 > (see (25)), (80) gives: 

dH H 2 T , 

— > A - 47rm 2 $ 2 (81) 

Consider once more the Hamiltonian constraint (59) which gives, using (79) 
to express k^k^ and (23) to express T 00 : 

\h 2 - 2 A - 87rm 2 $ 2 = WttU + 16tc{t 00 + pu 2 ) + a id a ij - R. (82) 

But it proved in [10], [26], that for the models under consideration, we always 
have: R < 0; (82) then gives: 

tt2 

A - 4vrm 2 $ 2 > 0; (83) 

3 

(81) then implies: 

so that H is non-decreasing. We also deduce from (83) since by (54) we have 
$ > 0, then $ > $(0) > 0, that: 

H 2 > 3A + 12vrm 2 ($(0)) 2 . (85) 

But by hypothesis, the r.h.s of (85) is strictly positive. So, since H is continu- 
ous, (85) implies: 

H < - V / 3A + 12vrm 2 ($(0)) 2 or H > ^3A + 127rm 2 ($(0)) 2 . (86) 

Also by hypothesis, H(0) < 0, then only the first inequality in (86) holds; 
moreover, (84) implies H > H(0) and (77) follows. □ 

We now introduce the Cauchy or initial value problem, taking into account 
(33), (35), (39); let the following quantities called initial data be given: 

g° = (gfj) a positive definite 3x3 constant matrix; 
k° = (k^) a symmetric 3x3 constant matrix; 
F° = (F^) an antisymmetric 3x3 constant matrix; 
u° = (u(°)' 1 ); E° = (E ' 1 ), constant vectors; 
$° > 0; U° > 0; p° > 0; e° > 0, real numbers; 

we look for g = ( gij ), k = (%), E = (E { ), F = (F^), u = (u*), $, U, p, e 
solutions of the evolution system such that: 

f g(0) = g°; k(0) = k°; E(0) = E°; F(0) = F°; «(0) = « (0) 

\ $(0) = $°; U(0) = U°; p(0) = p°; e(0) = e°. (87) 



Einstein-Maxwell-Massive Scalar Field System in 3+1 formulation 



15 



By Proposition 2.3, the constraint equations (59), (60), (61) ,(62) are satisfied 
if and only if the initial data satisfy these constraints we call initial constraints. 
In what follows, we consider that it is the case. 

We end this section by the useful notion of relative norm. Define the norm of 
a n x n matrix A by: 

|| A|| = sup { ^,i6K",^0}. 

If A\ and A 2 are two symmetric matrices with A% positive definite, define the 
norm of A 2 with respect to A\ by: 

IMMU^supj |gf, xgi",^o}. 

We have the following results proved in [17]: 
Lemma 2.5. 

\\A 2 \\ < H^lkPiH. (88) 
\\A 2 \\ Al < [tr(A^A 2 A^A 2 )]^. (89) 

Now by setting A\ = (g l:i ), A 2 = (ajj) a direct calculation gives: 

tr{A^A 2 A^A 2 ) = a i] a ir 

We then deduce at once from (88), (89), (90) that in these case: 

P 2 || < || A x || {a^aS- 

Next, let A = (a^) be a n x n matrix; set \A\ = sup-fla^l, i,j 
Then we have the following result, from [16]: 

Lemma 2.6. 

Let (u a ) = (u°,u l ) where vP and u l are link by (11). Then there exists a 
constant C > such that: 

<C\g\ 3 *; \F M \<(g lm F ol F 0m ^\g\l (92) 



(90) 

(91) 
= l,...,n}. 



u 
vP 



3 Local and global Existence of solutions 

We use an iterative scheme. 
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3.1 Construction of the iterated sequence 

We adopt the notations introduced in paragraph 2. We construct the sequence 
v n = (g n , K, E n , F n , u n , $ n , U n , p n , e n ), n G N as follows: 

• Set g = g°; k = k°; E = E°; F = F°; u = u®; $ = U Q = U°; 
Po = P°; e = e°. 

• If 9ni km F n , F n ^ 

Un, U n , p n , e n , are known: define T n a B, r njQ/3 by 
substituting g, F°\ F^, $ in the expressions (5), (6) of T a p, f a p by g n , 

Fn, F n , 

• Define v n+ i = (g n +i, k n +i, E n +i, F n+1 , u n+ i, $„+i, U n+1 , p n +i, e n+1 ) as so- 
lution of the linear ordinary differential equations (o.d.e) obtained by 
substituting g, k, E, F, u, $, U, p, e, T a/3 , T a/3 in the r.h.s of the evo- 
lution system (49) to (57), by g n , k n , E n , F n , u n , <f> n , U n , p n , 

It is very important to notice that, for every n the initial data for the linear 
o.d.e's are the same initial data g°, k°, E°, F°, u^°\ $°, U°, p° and e°. We 
obtain this way a sequence v n =(g n , k n , E n , F n , u n , $ n , U n , p n , e n ) defined in 
a maximal interval [0, T n [, T n > 0. 

3.2 Boundedness of the iterated sequence 
Proposition 3.1. 

There exists T > independent of n, such that the iterated sequence 

v n —(9n, k n , E n , F n , u n , $n ; U n , p n , e n ) is defined and uniformly bounded 

over [0, T[. 

Proof. 

Let AT G N, AT > 1, be an integer. Suppose that we have, for n < N — 1, the 
inequalities: 

9n - go\ < A\ ; (detg n )~ l < A 2 ; \k n - k \ < A 3 ; \E n - E \ < A A 
F n ~ F Q \ < A 5 ; \u n - u \ < A 6 ; \4> n - 0o| < A 7 
U n - U \ < A 8 ; \p n - pol < A 9 ; |e (n) - e (0 )| < A 10 (93) 

where A* > 0, i — 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 are given constants. 

We are going to prove that one can choose the constants A{ such that (93) 
still holds for n = N on [0,T[, T > 0, sufficiently small. Notice that the 
expression of (g%) = (g^ij)^ 1 contains (detgn)" 1 . 
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Integrating over [0,t], t > 0, the linear o.d.e satisfied by: g^, k N , E N , 
F N , u N , $at, U N , p N yields: 

f \9n - 9o\ < B t t ; \k N - ko\ < B 3 t ; \E N - E \ < B 4 t ; \F N - F \ < B 5 t 
I 107V - M < B 7 t ; \U N - U Q \ < B 8 t ; \p N - p \ < B 9 t (94) 

where Bi > 0, i — 1, 3, 4, 5, 7, 8, 9 are constants depending only on the 

We now study the cases of (detgN) -1 , un and e/v- 

The iterated equation satisfied by writes, using (49): 

9N,ij — — 2A;jv-i,ij- (95) 

Recall the formula: 

^[ln(detg N ] = g N d t g Ntij ; (96) 
on the other hand we have: 

^[ln(detg N )} = (detg N )~ l ^-(detg N ) = -(detg N )^(detg N )~ l . (97) 
(97) and (96) then give, using (95): 

— (detgN)- 1 = (2g N k N - hij )(detg N )~ 1 
an o.d.e in (deig^) -1 which integrate at once over [0, t], t > to give: 
{detgNY 1 = (detg )- 1 ^ [J 2(g%k N - 1 , ij )(s)ds) (98) 

Now (95) which is analogous to (49) shows that g^ and k^-i are the 
first and second fundamental forms of a space-like hypersurface; so 
g^fkN-i,ij = k % N _ li = £r(fcjv_i). We then deduce from (98) using (94), 
that: 

(detgN)- 1 < (detg )- 1 exp(dt) (99) 

where C\ > is a constant depending only on Ai and \g |, Hence, 
using (99), it appears that if we take in (93): A 2 > (detg )' 1 i.e 
(detg°)A 2 > 1, then given the continuity of t i — ?• exp(Cit), we will 
have for t sufficiently small (detg°)A 2 > exp(Cit) > 1. Then, there exits 
ti > such that, for < t < t\ we have, using (99): 

{detgNY 1 < A 2 - (100) 
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Next, in order to have for ujy and inequalities analogous to (94), 
we need to bound — which appears in the iterated equations in u n+ i, 

Pn 

e n+ i, built from equations (53) and (57) in u and e. Now following the 
definition of the iterated sequence (v n ) and given equation (56) in p, p n+1 
satisfies, for < n < N — 1 



Pn+l 



G r , 



where 



U n U n 



+ G 



t3 u° 



p n + Ane n g n . 



(101) 



(102) 



We have in the expression (102) of G n , using (92): 



u. 



u. 



< C\g n \i, and 



\u. 



1 u 

< 



II 



since u„ > 1 



Then, by (93), G n is bounded, i.e, 3 C > 0, \G n \ < C; hence, \p n +i\ < C and 
this implies: 

Pn+l > -C. (103) 

Integrating (103) over [0, t], t > yields p n+i > p° — Ct but since p° > 0, there 
exits t 2 > such that for t e [0, t 2 ] we have Ct < then 



(0 < t < h) (p n+1 > 



thus, < 4^ and the sequence 

Pn + l P 



Pr. 

is bounded. Hence un and ex also satisfy: 

|«at - u Q \ < B 6 t; \e N -e \<B w t 



(104) 



where as in (104) B 6 , B w are constant depending only on the A\. We then 
conclude that if T > is such that: 



T<inf(t 1 ,t 2 ), B{T<A U i = 1,2,..., 10, 

then, by (94) and (104) vn also satisfies (93) on [0, T[. Hence, the iterated 
sequence (v n ) is uniformly bounded over [0,T[. □ 
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3.3 Local existence of solutions 
Theorem 3.2. 

The initial value problem for the Einstein- Maxwell- Scalar Field system has a 
unique local solution. 

Proof. 

We prove that the iterated sequence (v n ) converges uniformly on each bounded 
interval [0,5] C [0, T], S > 0. For this purpose, we study the difference 
v n +i — v n . But given the evolution equations (54) and (55) in <3> and U, we will 
deal with the difference: 

^2U n+1 - y/U n - ~ 7m ^ Wn 

We then need to show first of all that the sequence 

is uniformly bounded. 

• By (55), the iterated equation providing U n +\ writes: 



U n+1 = 2H n U n - mH n ^2~U n . (105) 

But by Proposition 3.1, there exits a constant C > such that we have 
over [0, T[: 

\2H n U n - m 2 $ n Jw n \ < C 



(105) then gives: 

dU n+1 



> -C 



dt 

and integrating over [0,i], < t < T yields: 

U n+1 >U°- Ct. 

Recall that U° > 0. Then taking t sufficiently small such that Ct < ^ 
we have U n+ i > Then 

1 1 



which shows that y v ^r ) is uniformly bounded over [0, T[, T > small 
enough. 
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Since (— ) is also bounded, taking the difference between two consecutive 
iterated equations, we deduce from the evolution system, using 
^n(O) = v°, Vn, that there exits a constant C2 > such that: 

IflW-i ~~ 9n\ + \k n +i — k n \ + \E n+ i — E n \ + |F n+ i — F n \ + \u n+ i — u n \ + 
|$ n+ i - $„| + \U n+ x - U n \ + \p n+ i - p n \ + \e n+ i - e n \ < 

C2 J (\g n — 9n-i\ + \k n — k n -i\ + \E n — E n -i\ + \F n — F n _i| + 

\u n - U n -l\ + |$ n - $ n _l| + |f/„ - C/ n _l| + \p n ~ p n -l \ + 

|e„ -e n _i\)(s)ds. (106) 



So, if we set: 

= Ifl'n+l — 9n\ + |^n+l ~ fc n | + \F n+1 — F n \ + \E n+1 — E n \ + \u n+1 — U n \ + 

- $ n | + \U n+1 - U n \ + \p n+ i - p n \ + \e n+1 - e n \ (107) 
then (106) shows that we have: 

On(t) < ^2 /* On_i (s)ds. (108) 
JO 

(108) gives, by an immediate induction on n: 

«n(*) < ll^llocf 2 ^- (109) 

[n — 2)1 

But since the series J2 is convergent we have necessarily a n (t) — > 

n=0 

as n — > +oo. Definition (107) of a n then shows that each of the 
sequences (g n ), (k n ), (E n ), (F n ), (u n ), ($ n ), (U n ), (p n ), (e n ) converges 
uniformly on each bounded interval [0,5], < 5 < T and that their 
respective limits denoted: g, k, E, F, u, U , p and e are continuous 
function of t. From the iterated equations, it appears immediately that 
there exits a constant C3 > such that: 



dv n+ i dv n 



dt dt 



<C 3 \v n+1 -v n \ (110) 



where we set, for v = (v l ), \v\ = J2 I v% I- The convergence of (v n ) then 
implies, given (110), the convergence of 

dv n 

In 

and hence, that each of the sequences (g n ), (k n ), (E n ), (F n ), (w„), ($n), 
(U n ), (p n ), (e„) converges uniformly on each interval [0, 8], < S < T. 
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Consequently the limit functions g, k, E, F, u, $, U, p and e are of class 
C 1 and v := (g,k, E, F,u,$,U, p,e) is a local solution of the coupled 
Einstein-Maxwell-Massive Scalar Field system. 

Now to prove the uniqueness of the solution, suppose V\ and v 2 are two 
solutions of the Cauchy problem , with the same initial data. Then, 
defining a{t) the same way as a n (see (107)) for the difference \v\ — v 2 \, 
leads, using the evolution system to: 

a(t) < / a(s)ds 
Jo 

which gives by Gronwall Lemma a = 0; hence v\ = vi and uniqueness 
follows. 

□ 



3.4 Global existence of solutions 

We prove: 
Theorem 3.3. 

Let A > — 47rm 2 ($ ) 2 be given and suppose H° := (g )*- 7 ^ < 0. Then the 
initial value problem for the Einstein-Maxwell-Massive Scalar Field system has 
a unique global solution defined all over the interval [0, +oo[. 

Proof. 

Following the standard theory on the first order differential systems, it will be 
enough if we could prove, given the evolution system (49) to (57) that, if each 
of the functions: \g\, \k\, \E\, \F\, |$|, \u l \, |p|, |e|, \Rij\, \R\, (detg) -1 , \-\, is 
uniformly bounded over every bounded interval [0,T*[, where T* < +oo. 
Notice that the hypothesis of Theorem 2.4 are satisfied; so (77) applies, i.e H 
is bounded. 

• We deduce at once from (83) and using (77) that $ is bounded. 

• Since the l.h.s of (82) is bounded and since —R > 0, r 00 > 0, a^a^ > 0, 
U > 0, p > 0; we deduce that U and pu\ = p(u ) 2 are bounded; but 
p( u o) 2 > P > 0, since u° > 1; then we deduce that p is bounded. 

• (80) gives, using (79) to express cr y f7y, and since: U > 0, g^UiUj > 0, 
pu\ > 0, Too > 0: 

ATT 

> kak ij - A - 4vrm 2 $ 2 . (Ill) 

Then integrating (111) over [0, t], < t < T*, we have, since H and $ 
are bounded: 

/ k ij k lJ {s)ds < +oo (112) 
Jo 
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Next we have, integrating (49) over [0,t], t G [0, T*\. 

\9{t)\< 1/1 + 2 f\k{s)\ds (113) 
Jo 

but setting Ai = (g lj ), A 2 = (hj), (91) gives: 

||A:|| < WgMhj*?*)* (H4) 

and we deduce from (113) and (114) that: 

\\g(t)\\ < \\g°\\+2 [ t \\g(s)\\(k ij k i ^(s)ds. 
Jo 

Hence, by Gronwall Lemma, there exists a constant C > such that: 

\\g(t)\\ < C\\g°\\ex V (C f\k l3 k l ^(s))dt (115) 



but we deduce from (115) applying Schwarz inequality, using (112) and 
since T* < +oo, that \\g\\ and hence \g\ is bounded. 

By (58), Rij expresses in terms of 7^ given itself by (19), which involves 
g l i; so we need to control (detg) -1 . We use once more the formula: 

Then, using the evolution equation (49), we obtain: 

j t [Hdetg)] = -2H. (116) 

Since H is bounded, we obtain, by integrating (116): 

-C < \n(detg) < C 
where C > is a constant. Hence: 

e~ c < detg < e c 

which shows that, both detg and (detg)^ 1 are bounded. Then, by (58), 
Rij is bounded and R = g^Rij is bounded. 

Now consider the Hamiltonian constraint (59) which gives, since Too > 0, 
pu\ > 0, Too > (see(23)) 

R + H 2 -2A> k {j k ij . 

Then, since R and H 2 are bounded, kijk^ which is positive is bounded. 
Then, since \\g\\ is bounded, by (114) , \\k\\ and hence \k\, is bounded. 
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Deduce from (82) whose l.h.s is bounded and using U > 0, r 0Q > 0, 
pul > 0, <J % i<Jij > 0, —R > 0, that r 00 is bounded. But by (25): 



which implies: 



too = \gv&E? + 



< ^g ij E i E j < roo ; < -^F^ < r 00 



(117) 



:il8) 



and gijE l Ei is bounded. Hence, using (92) and since \g\ is bounded, E l 
is bounded. 

The constraint (62) shows, since u° > 1 and E % is bounded, that e is 
bounded. 

Integrating the equation (52) in F^ over [0, t), t < T* < +oo shows, since 
\g\ and \E\ are bounded, that \F\ is bounded. 

It remains the cases of - and u % . 

p 

Expression (32) of p gives, using the notations (48): 



u 



exp 



\ (-H + Clj-^ds . (119) 
Jo J u u J 



p°u°(0) 

Now as we already indicated, equation (8) gives for (3 = 



u 



U U 



4neg, 



u l Ei \ 1 



13 



P 



We deduce from (119) and (120) that: 

ii° < G(t)u° 

where: 

|%wV| \^ireg i: jE j \ \u 

(*[t) — , nNn H 



0>|2 



exp 



p°u°(0) ' u° 
Integrating (121) over [O.t], t <T* < +oo gives: 

u°(t) < M°(0)exp [ f G{s)ds 



(-H + Cl- )(s)ds 



(120) 
(121) 
(122) 

(123) 



u 



then, using (92) to bound — and since k, e, g, E, H are bounded, by 

(123), there exists a constant C(T*) > such that u° < C(T*)u°(0). 

Hence by (119), - is bounded and writing u % = %u° shows that u % is 
p u 

bounded. This completes the proof of theorem 3.3 



□ 
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4 Asymptotic behaviour 



We consider the global solution over [0, +oo[ and we investigate the asymptotic 
behaviour of the different elements at late times. We introduce the following 
quantity which plays a key role: 



Q = H 2 - 247rTo - 3A 
At late times, we have the following asymptotic behaviour: 
Theorem 4.1. 



where 



A + 47rm 2 ($°) 



0\2 



(124) 



Q 




0(e 7 ) 


(125) 


0-ij(T ij 


— 


0{e 7 ) 


(126) 


\R\ 


— 


0(e 7 ) 


(127) 


T~00 


= 




(128) 


p 


= 


0(e 2lt ) 


(129) 


($)2 


= 


0(e 7 ) 


(130) 








(131) 


E*E t 




0(e- 27 *) 


(132) 


Ht)\\ 


< 


Ce- 7t |k(t)|| 


(133) 






L > 


(134) 


Too 




m 2 

Lj 

2 


(135) 


H 




-(3C )5 + O(e- 25t ) 


(136) 


~ 25 %\ 


< 


C 


(137) 


e 2St g ij \ 


< 


c 


(138) 


9ij(t) 




e 28 \G l3 + G(e- 7 *)) 


(139) 


9 ij (t) 






(140) 


\k--\ 

\ hj %3 \ 


< 


Ce 25t 


(141) 


\E l \ 


< 




(142) 


\F--\ 


< 


Q e (2S+v)t 


(143) 


\F ij \ 


< 




(144) 


e 


< 




(145) 



7 



v = 35 — 7; 



C = A + 4vrm 2 L; 5 



[-(A + 47rm 2 L)] 5; 



(146) 



Gjj a symmetric positive definite constant matrix; (G y ) = (Gy)" 
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Proof. 

Notice that by (23) and (35), we have: m 2 $ 2 = 2T 00 - 2U. Expression (124) 
of Q then shows, using (82) that Q > 0. Let us point out first of all that, the 
quantity Q defined by (124) plays a key role in GR, and in the presence of the 
massive scalar field, it stands for the quantities S in [13], Z in [20], S in [11], S 
in [7] and reduces to H 2 ± 3A in [26] which deals with the case of zero scalar 
field. 

• We have, using the expression (23) of T 00 and (35): 

T 00 = U+^m 2 ® 2 (147) 

then, the evolution equations (54) and (55) in $ and U give: 

T o = 2HU. (148) 

Expression (124) of Q then gives, using (148): 

Q = 2H(H - 2AnU) 

then, using equation (63) in H, in which we use the Hamiltonian con- 
straint (59) to express R + H 2 and (79) to express kijk^ and since by 
(23) and (35): Ang^T^ = 2AttU - 12T 00 , we obtain 



Q = -H 
^ 3 



iJ 2 -247rToo-3A+3(87rroo+47rpM 2 +47r^' J M i M j +a ii a ii ) . (149) 



But since 87TToo + ^npu^ + Ang^UiUj + cr^cr^ > 0, H < 0, and given the 
definition (124) of Q, (149) gives: 

Q < \HQ. (150) 
Integrating (150) over [0, t), t>0 yields: 

< Q < Q exp[f -tfds]; 

Jo 3 

and (77) gives 

< Q < Qo exp ( - 2tJ 1 (A + 47rm 2 ($ ) 2 ; 



3 



and (125) follows. 



Using (82), the results (126), (127), (128), (129) and (130) (since ($) 5 
2U) are direct consequences of (125). 
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(131) and (132) are consequences of (128), using (118). 
Setting A\ = (g^) and A 2 = (cr^), (91) gives: 

Ht)\\ < 

(133) then follows from (126). 

The evolution equations (54) and (55) in $ and U give: 



(151) 



m $$ + U 



d 
dt 



m 



2HU. 



(152) 



But since H < and U > 0, (152) implies: |[^$ 2 + U] < 0; we then 
deduce, using U > that: 



^$2 < + £/ < ^($0)2 + ^0 

2 ~ 2 - 2 v ; 



(153) 



Hence $ is bounded. But the evolution equation (54) in $ shows that 



$ > 0; then $ > $ u > and since 



2$$ > 0, $ J is an 



increasing function. $ 2 being positive, increasing and bounded has a 
strictly positive limit, i.e, there exits L > such that 



with: 



(154) 
(155) 



and we have (134). 

(135) follows from (147), (130) and (134). 

To prove (136) which is one for the main results, since by (125) Q 
its expression (124) shows, using (135) that: 

H 2 -3A — > 12vrm 2 L. 

Hence: 

H 2 -(3A+12nm 2 L) = #-(3A+127rm 2 L)5 #+(3A+127rm 2 L)^ 

But by (77), H < 0; so: 

H - (3A + 12nm 2 L)^ < -(3A + 12nm 2 L)^ < 0. 
We then deduce that: 



(156) 



H 



-(3C 



(157) 
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where: 

C = A + 4rrm 2 L > 0. (158) 
But since by (84) H is an increasing function (157) implies: 

H<-(3C )^. (159) 
Now (155): — 47rm 2 $ 2 > —47rm 2 L; (81) then implies: 

^>—-C (160) 
dt ~ 3 K ' 

write: 

^ - C = \{H 2 - 3C ) = \i-H + (3C )i][-# - (3C )*] 
in which using (159) we have: 

-# + (3C )3 > (3C )* + (3C Q )z = 2(3C )*; 
(159) also implies: —H — (3C )^ > 0. We then deduce from (160) 

d E>26[-H-(3Co)^} (161) 

where: 

5=l(3C )l (162) 

Write (161) in the form: 

j t (H + (3C )3) + 25[tf + (3C )5] > 0. (163) 

Multiply (163) by e 2St > and integrate over [0, t] to obtain: 

e 25 '[# + (3C )5] > H(0) + (3C )5. (164) 

Now multiply (164) by -e~ 25t < 0, use once more (159) which gives 
H + (3C )5 < to obtain: 

|#+(3C )^| < |iJ(0) + (3Co)^|e- 25t 

and (136) follows with , see (158) and (162): 5 = [|(A + 47rm 2 L)]i 



JV. Noutchegueme and A. Nangue 
To prove (137), set = e~ 25t gij] then we have, using equation (49) in 

ga- 
rth ■ ■ 

^l = ^26h lJ -2k l ,e- 25t . (165) 

Now, using (78) to express ky, the result (136), the expression of C and 
S in (146) we deduce from (165) that: 

rJh ■ ■ 

^ = 0(e~ 2H )h tJ - 2e~ 2 %, (166) 
Integrating (166) over [0, t], t > 0, and taking the norm yields: 

\\h(t)\\ < \\h(0)\\ +C [\e- 25s \\h(s)\\ +e- 2Ss \\a(s)\\)ds. (167) 
Jo 

where C > is a constant. Notice that W = e 2St g^. Now setting in 
formula (91): A x = (h ij ), A 2 = yields: 



W(s)\\ < WKsMa^ 



I 2 



where o% = h} l h? k o lk = e 45t a ij . Hence, ||o"(s)|| < \\h(s)\\e 25t (a^ a^- 
(126) then gives 

IkOOII < IIM^)l|e 25s e-^. 
We then deduce from (167): 

\\h(t)\\ < ||/»(0)||+C [\e- 2Ss \\h(s)\\+e^ s \\h(s)\\)ds 
Jo 

then, since 7 < 5 [see (146)]: 

\\h(t)\\<\\h(0)\\+C [\-^\\h(s)\\ds. 
Jo 

By Gronwall Lemma, this gives: 

< ||/»(0)||exp[/* C7e"^-] < C7, 
Jo 

where C > is a constant, and (137) follows. 

To obtain (138) set this time L u = e 2St g^ and proceed as for (136) and 
obtain (138). 

To prove (139) which is one of the main results, first use 

e - 25t aij = 0{e^ 1 ) (168) 
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which is a direct consequence of (151), (126) and (137). Recall that 
setting hij = e~ 26t gij led to (166). We deduce from (166), using (168), 
7 < 5, and since h is bounded: 

^ = 0(e-T*). (169) 
at 

(169) shows that has an exponential fall of at late times and by the 
mean value theorem, hy has a limit we denote G^ as t — > +oo. Then 
we can write: 

fn j (jt) = G ij + 0{e^ t ) (170) 

where, given the properties of hij, Gij is a symmetric, positive definite 

in- 



constant 3x3 matrix. (139) follows from (170) since = e 2St hi 



• (140) is a direct consequence of (139) since (g^) = (g^) 1 . 

• To prove (141), use (78) which implies, since H is bounded 

\\k\\<C(\\g\\ + \\a\\). (171) 
(141) then follows from (171), using (139) and (168). 

1 3 

• To obtain (142), use (92) which gives \E l \ < [E 3 Ejp\g\^ and conclude 
by applying (132) and (139). 

• To obtain (143), integrate equation (52) in and use (139) and (142). 

• To obtain (144), use F ij = g il g jk F kl , (140) and (143). 

• Finally to obtain (145), use the constraint equation (62), (142) and 
u° > 1. This completes the proof of Theorem 4.1. 

□ 

Remark 4.2. 

The result (139) shows an exponential growth of the metric tensor g at late 
times, this result, but also (136) confirm mathematically the accelerated ex- 
pansion of the universe as observed in Astrophysics. 

5 Geodesic Completeness 

We prove: 
Theorem 5.1. 

The space-time which exists globally is future geodesically complete. 
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Proof. 

We use the fact that, the geodesies equations for the metric (1) imply that, 
along the geodesies whose affme parameter is denoted by s, the variables t, 
u°, u % satisfy a first order differential system containing between others, the 
equation: 

Ts = "° < 172 > 

The space-time will be future geodesically complete if we prove that the affme 
parameter s also goes to infinity. Then using (11), the notations (48) and 
(172), it will enough if we could prove that: 

^ = (l + ft) -«V)-bC>0 (173) 



where C > is a constant, since one could then deduce at once from (173), 
integrating, that s > Ct + D where D is a constant; hence s — > +oo as 
t — > +oo. Our goal will then be to prove, that (1 + gijU l vP) or finally gijU l u J , 
is bounded by a strictly positive constant. For this purpose we use equation 
(53) in u % . It shows to be useful considering Ui = gijU 1 , rather than u % . Differ- 
entiating this relation, we have: 

dm du J Agij 

— — = Qij— -, — ha J -r (174) 
dt y j dt dt ( J 

then, using equation (53) in u % and equation (49) in g^, we deduce from (174), 
the equation: 

du l ,• u l u k , e : e „ u l . , 

-e = -wA-sr + »(- V + V V- < 175) 

In order to bound (gijtfvP) we set up a differential equation for this quantity. 
First notice that equation (49) gives, using g^gu = 5j : 

from where we deduce: 

A direct calculation using (175) shows that, in (176) we have: 

,-„■ dui e . e „ u l u k ,■ u^u k Ui , 
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But given the antisymmetry of Fi k , we have F[ k uu = 0. Next using the 
expression (19) of 7' - we obtain for the last term in (177): 



Y jk U j U k Ui 



-lgjiC l kp u j u k u p + lg kl C l pj u j u k u p + \c) k u 3 u k 



U: 



But since C l j k = — C k j the last term in (178) vanishes. For the same reason 



(178) 



1 



1 



1 



-29iiC l kp u j u k u p + -g kl C l pj u j u k u p = -(C l pk + C l kp )g 3l u 3 u k u p = 0. 

Consequently, the r.h.s of (177) reduces to its first term and (176) gives: 

d ■• ■• e 

— (g^UiUj) = 2k v UiUj - 8ir-UiE\ 

Ctv P 

But by (78), k i3 = fg i3 + <7«; (179) then gives, using (136) and (146): 
d 



(179) 



dt (9«u iUj ) 



25 + ©(e-^J^UiMj + 2<r %3 u i u j - 8n- Ui E l . 



;iso) 



Now we have a 13 = g tk g 3l o~ k i, so by (168) and (140): e^+^V 1 - 7 is bounded. 
This implies, since the matrix G lJ is positive definite and constant, that there 
exits a constant C > such that: 



o~ i3 UiUj < Ce-^+^G^muj. 



Now by (140), there exists a constant C > such that: 

G ij u iUj < Ce 2St g ij Ui u r 
Now deduce from (32) and (38) that: 

e _ e° 

we then obtain from (180), using (181), (182) and (183): 

~(g ij u iUj ) < (-25 + Ce- 2<5 VWj + Ce^g* muj + C\ Ui E l \. 

Ctv 

Now set: 

then we have, using (184): 
dW 



W = e 25t g ij u iUj 



(181) 
(182) 

(183) 

(184) 
(185) 



d 



25e 25t g ij u iUj + e 2dt ~{g tJ u iUj ) 



dt dt 

< Cg'u.uj + Ce^-^g^UiUj + Ce 2St \ Ul E l \. 



;i86) 
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Now since g is a scalar product: \ii{E l \ < {u l Ui)^{E l Ei)2 . so, by (132) we have 
\uiE l \ < C(g lj u i u i )^e~' rt ; (186) then gives: 



dW 
dt 



< Ce~ 25t W + Ce-TW + Ce^^s 



^187) gives, since < 7 < 5: 

dW 



dt 



< Ce'^W + Ce 5t wl 



But it is well known that by (188) we have: 

W{t) < z(t) 

where 

I = Ce-^z + Ce bt z\ 
z(0) = W(0). 

But (190) is a Bernoulli equation whose solution is: 



(187) 

(188) 

(189) 
(190) 



z{t) =exp(-2^*a(s)rfs) (W(0))3 + J* b(s)(exp J\(r)dr)ds 2 (191) 



where: 



a(t) 



(192) 



One deduces easily from (191), (192) that: 

z(t) < Ce 2St , 

where C > is a constant. Then using expression (185) of W and (189), we 
obtain: 

gijU l u J < C. 

This completes the proof of Theorem 5.1. □ 



6 Energy conditions 

In this section we prove that the global solution satisfies the weak and the 
dominant energy conditions and, under some hypothesis, the strong energy 
condition. Recall that a viable physical theory is supposed to fulfill at least 
one of the energy conditions (Hawking[ ]). In fact notice that considering the 
stress-energy-matter tensor of the Einstein equations (2), and keeping the (~) 
to avo'id any confusion, the quantity 
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represents physically, the energy density of the charged particle, measured by 
an observer whose velocity is V a and so must be non-negative, V a being a 
future pointing time-like vector, see [27]. 

We recall below the three types of energy conditions: let (V a ) and (W a ) be 
any two future pointing time-like vectors. The solution is said to satisfy: 

1) the weak energy condition if: 

{f a p + r a/3 + pu a up)V a V? > 0. (193) 

2) the strong energy condition if: 

R a pV a V f) > 0. (194) 

3) the dominant energy condition if: 

(T aP + r afi + pu a u fi )V a W p > 0. (195) 

Obviously, (195) implies (193), just setting: W a = V a . 
We begin by proving: 

Proposition 6.1. 

Let V a and W a be two future pointing time-like or null vectors. Then 

V a W a < (196) 

Proof. 

Since V a W a = g a j3V a W 13 and given the definition (1) of g, (196) is equivalent 
to: 

- V°W° + gijVW < 0. (197) 
Now, V a , W a being future pointing time-like or null we have: 

V a V a < 0; V°> 0; W a W a < 0; W > 

,or, equivalently: 

< gijV^ < (V ) 2 ; V° > 0; < gijWW < (W ) 2 ; W° > 0; 

hence: 
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which gives: 

< (g lj V i V j )^g ij W i W j ) 1 ^ < V°W°. (198) 
But since g is a scalar product, we have: 

\gi^rW'\ < {g l3 V l V^{g tj W l W^ (199) 

(197) then follows from (198) and (199). □ 

Next we prove this important result for the Maxwell tensor r a p defined by 

(6): 

Proposition 6.2. 

For any two future pointing time-like vectors {V a ), (W a ), we have: 

T a pV a W a > 0. (200) 

Proof. 

It will be enough if we could prove (200) by choosing any suitable frame. Let 
us consider the frame of the four vectors: 

I = (l a ); n = (n a ); x = (x a ); y = (y a ), satisfying the following properties: 

lj a = n a n a = l a x a = l a y a = n a x a = n a y a = 0. (201) 

Now inspired for instance by the case where the electromagnetic fields F a p 
derives from a potential vector, the antisymmetric 2-form F a p can be written 
in one of the two following general forms: 

~ A B 

F a p = —{l a np - lpn a ) + -w(x<xyp ~ %pya) (202) 

or: 

~ C 

F a j3 = ~ lpX a ) (203) 

where A, B, C are constants. In fact, since M = K x G, with G is a simply 
connected Lie group, by Poincare Lemma, F a p is an exact form. Then, there 
exits a potential vector A a over M such that: F a p = V ' a Ap — V ' pA a then 
formula (202) generalizes the form given by the development of the above 
expression of F a p by applying (21) in the case of the frame (e a ) whereas (203) 
corresponds to the case of the natural frame (d a ). 

• Next, it shows useful to choose the constants A, B, C by assuming that 
we have in addition: 

l Q n a = -l; x a x a = y a y a = 1; x a y a = (204) 
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Now consider the Maxwell tensor (6) i.e 



Tap = -\g a pF x »F Xll + F a \Fp x (205) 



1°) Consider the form (202). 

A direct calculation using (201), (204) gives: 

~ ~ B 2 - A 2 ~ ~ A 2 B 2 

F Xti F Xll = ; FaxFp A = — (l a np + njp) + — (x a Xp + y a y p ) 

(206) 

so that, in this case, (205) and (206) give: 



1 



(A 2 — B 2 \~\ 

A 2 (l a rip + njp) + B 2 (x a xp + y a yp) +g a py J . (207) 



Now if we express the vectors V = {V a ), W = (W a ) in the frame 
(l,n,x,y) by: 

( V a = Ml a + Nn a + Px a + Qy a (208) 
[W a = MT + N'n a + P'x a + Q'y a (209) 

and if we set: 

Kp = -l a np - njp + x a x p + y a y p (210) 
then a direct calculation, using (201) and (204) gives: 

KpV a W p = g aP V a W p = (-NM' - MN' + PP' + QQ'). 

Hence g a/ 3 = h a p. So, we can express g a p in (207) by (210) and this 
gives: 

A 2 + B 2 ^ . . 

Tap = g (211) 

where 

T a p = l a np + njp + x a xp + y a yp (212) 
• Now a direct calculation using (208), (209), (212), (201) and (204) gives: 

7 a pV a W f) = NM' + MN' + PP' + QQ. (213) 

But if V — iy a ) and W = (W a ) are future pointing time-like vectors, in (208) 
and (209), we add: 

V a V a < 0; W a W a < 0; M > 0; M' > 0. (214) 

Now (214) writes, using (208), (209), (201) and (204): 

-2MN + P 2 + Q 2 < 0; -2M'N' + P' 2 + Q' 2 < 0; M > 0; M' > 0; 
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so we have: P 2 + Q 2 < 2MN; P' 2 + Q' 2 < 2M'N'\ M > 0; M' > 0. But this 
implies since M > 0; M' > 0: 

p2 i o 2 P' 2 + O' 2 
V > 0: \ ' > (I: MJV > ^ - : .\/'.V > ^ " . (215) 

Then, since M > 0; M' > 0, (215) gives: 

iVM' > — (P 2 + Q 2 ); N'M > —(P' 2 + Q' 2 ); (216) 
~ 2M V - 2M' 

So if we consider (213) in which V = (V a ) and W = (V a ) are future pointing, 
we deduce from (216) that: 



2MM' 



M' 2 (P 2 + Q 2 ) + M 2 {P' 2 + Q' 2 ) + 2MM'{PP' + QQ') 

(217) 

Considering the term in the square bracket in the r.h.s of (217) as a quadratic 
polynomial in M, its discriminant is: 

A = M' 2 [(PP' + QQ') 2 - (P 2 + Q 2 )(P' 2 + Q' 2 )]. 

But by the properties of the usual scalar product in IR 2 : 

(PP' + QQ') 2 < (P 2 + Q 2 )(P' 2 + Q' 2 ) 

then A < and the r.h.s of (217) remains positive, and so is the l.h.s. Then 
in this case, (200) follows from (211). 

2°) Consider the form (203). 
A direct calculation using (201) and (204) gives this time: 

~ ~ ~ ~ C 2 

F Xfl F\n = 0; F a \Fp x = —l a lp. 

Then (205) gives 

C 2 

T aP = —IJp (218) 

so if V = {V a ) and W = (V a ) are two future pointing time-like vectors, using 
the decomposition (208) and (209), we obtain from (201) and (204): 

TapV-Wp = ~(l a V a )(lpWe) = —(-AO(-iV') = °^NN'. (219) 

But (215) which holds since V and W are future pointing, gives N > and 
N' > 0. Hence, by (219) r a pV a W^ > 0. This completes the proof of Proposi- 
tion 6.2. □ 
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Theorem 6.3. 

The global solution of the coupled Einstein-Maxwell-Scalar Field satisfies: 
1°) the weak and the dominant energy conditions. 

2°) the strong energy condition if A > . 
Proof. 

1°) we first prove that the solution satisfies the dominant energy condition 
(195). Let V = (V a ) and W = (W a ) be two future pointing time-like 
vectors. By (200) we have 

T a pV a W f> > 0. (220) 

Next we have, since u = (u a ) is a time-like future pointing vector and 
using (196): 

{pu a u p ){V a W p ) = p{u a V a ){u p WP) > 0; 

so: 

{pu a u p ){V a W p ) > 0. (221) 
Now the expression (23) of T a p gives: 

T a pV a W fi = f 00 V°W° + TijVW 

= ^($ 2 + m 2 $ 2 )V°W° + ^($ 2 - m 2 $ 2 )^VW J ' 

then: 

T ap V a W p = ^ 2 {V°W° ' + g^W^ + ^m 2 ® 2 ^ ^ - gij V l W j ). (222) 

But by (196) which is equivalent to (197), the last term in the r.h.s of 

(222) is positive. 

Next, since g is a scalar product and V, W are future pointing vectors, 
we deduce from (198) and (199), that 

IgijVW^ < V°W°; 

then 

g ij V i W j > -V°W°. (223) 

(223) then shows that the first term in the r.h.s of (222) is also positive. 
Consequently, 

T a pV a WP > (224) 

(195) then follows from (220), (221) and (224). Hence the dominant 
energy condition (195) is satisfied. 
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• Setting in (195), W = V, we have (223). Hence the weak energy 
condition (193) is satisfied. 

2°) We now prove the strong energy condition (194). 

• Let V = (V a ) a future pointing time-like vector. We deduce from 
the Einstein equations (2) that: 

RapV a V^ = {~R - A)V a V a + 8vr(T a/3 + f Q/3 + pu a u p )V a V p . (225) 

Since (193) is satisfied, the second term in the r.h.s of (225) is positive. 
In the first we have : 

R = g a0 R a /3 = 9°°Roo + 9 l3 R l3 = -Roo + 9 l3 % r (226) 

Recall the classical formula linking R^ and Rif 

Rij = Rij dfkij + Hkij Iknky 

Contracting by g l i yields: 

g ij Rij = R- g^dthj + H 2 - 2hjk ij (227) 

write: g lj d t kij = d t (g l3 kij) — kijd t g lj = d t H — 2kijk' lj , (227) then gives: 

g ij Rij = R-d t H + H 2 . (228) 

Now setting in the Einstein equations a = (3 = yields: 

~ R 

^oo = -- + A + 8vr(Too + r o + pMo)- (229) 
(226) gives, using (227) using (229): 

R=^-A- 87r(T 00 + fbo + pul) + R-d t H + H 2 . 

From where we deduce, using T 00 + r 00 + pu^ > 0, R < 0, —d t H < 0, 
(given by (84)): 

R 

— - A < -2A + H 2 . (230) 
But by (77) we have: H 2 < (#(0)) 2 . Hence (230) gives: 

- - A < -2A + (H(0)) 2 (231) 
2 
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but by hypothesis: -2A + (H(0)) 2 < 0; hence: 

f-A<0; (232) 

since V a V a < 0, (232) implies that the first term in the r.h.s of (225) is 
positive; we conclude that; we have: R a pV a V^ > 0. This completes the 
proof of Theorem 6.3. 

□ 



Concluding Remarks 

In our future investigations, we will take into account the aspect "distribution" , 
of the charged particles. For this purpose we will couple the Vlasov (resp. 
Boltzmann) equation in the collisionless(resp.collisional) case. 
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7 Appendices 

A\. Proof of formula (22) 

We use the Codazzi equations which write: 

Rxi,jin x = -ViKj + Vjk u . (233) 
But n = (1,0,0,0) so (233) gives: 

Roi,ji = -Vihj + Vjk u . (234) 
Now the curvature tensor, on (M,g) writes, see [2], p. 240: 

R\,^ = - M^o) + H^o) - ill^a) - C^lL- (235) 

In particular, taking in (235): A = I; a = j; (3 = i\ fi = 0, we obtain, using 
(21), (12), (13), (18) and (19): 

R l ^ = -j t {l\ j )-Vik l j . (236) 

Now (234) gives, using the symmetry properties of fl A ft ^: 

V% • V,/.-!. (237) 
Equalize the two values of R l j ji0 provided by (236), (237) to obtain (21). 

A2. Proof of Lemma 2.1 

1°) Proof of (63) . 

First deduce from the evolution (49) in g^\ using g d gji = 5* that: 

^ = 2k lJ . (238) 
dt v ; 

Now since H = g^k^, (238) gives: 

f^fw- < 239 > 

A direct calculation using the evolution equation (50) in k^ and the 
relation g al3 f a p = which implies g^Tij = t 00 , gives 

dk- ■ 

gV ^df = R+H *- 2k V kK '+W J \T ij + P u i u j )-l2ix(T 00 +pul)-9>TXT m -'iK. 

(240) 

then (63) follows from (238) and (240) 
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2°) Proof of (64) 

We have ^ = 2H^§; then use (22) to obtain (64). 



3°) Proof of (65) 
We have 



d(kijk tj ) dk lj jijdky 



M k *lU+ k M- < 241) 
We also have: k^ = g ll g^ m ki m \ so, we have, using (238) 

d] ^ = Ak a kl+g a g^ d ^ (242) 
(242) gives by a direct calculation, using the evolution equation (50) for 

klm- 

— ■ - R ij + Hk ij + 2kjk il -87r(T ij + r ij + pu i u j ) 



dt 



+ 4ng 



T 00 -pul + g lm (T lm + p Ul u m ) -Ag ij . (243) 



(65) then follows from (241), (243), using once more the evolution equa- 
tion (50) for kij to express the last term in (241). 

4°) Proof of (66) and (67) 

We use the conservation laws: 

V Q (T Q/3 + f af} + pu a u p ) = (244) 

(244) writes, using the formulae (21) 

e a {f aP +T aP +pu a u p )+^x{T xp +T Xp +pu^ = 0. 

(245) 

It shows useful to write (243) in the form: 

e (T ^ + r ' 3 + puV) + e,(^ + t 1 ? + pu l u?) + ^ (T ? + + puV) 

(246) 

Then: set in (246) (3 = and use (19) to obtain (66). 
in (246) (3 = j and use (19) to obtain (67). 

5°) Proof of (68) 
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We have R = g^Rij, where Rij given by (58), then: 

ai 

So we have, using (238) 



dR _ dg^_ ij dRg 
dt ~ ij dt +9 dt 



2R ij k ij + g ij ^-. (247) 



Now we verifies by a direct calculation, using the formula (21) and the expres- 
sion (58) of Rij that: 



dRij _ (d% 



(248)gives, using formula (22) and since H = g^kij depends only on t: 

<7<^ = -2V,V,* fl (249) 
(68) then follows from (247) and (249). 

A 3 . Proof of (70) 

It is easily seen, using the definition of Aj in Lemma 2.2, that 

Ai = 9j iB j (250) 

where: 

Bi = Vik ij - 8Tr(T 0j + r 0j + pu°u j ). (251) 
We then have, differentiating (250) and using equation (49) in g^: 

dAi , * dB j 

- = -2^ +J ,— . (252) 

Now we have by (251): 

dRi d d 

—jr = |r(W) - StAt* + r« + p«V) (253) 
dt dt dt 

But Vik 1 ^ = dik^ + Yuk^ + Jiik l1 ; then we have: 

Now use (243) to express (22) to express and obtain: 

d{ yjfl = y . R ij + HV . k H + 2 (V \k ll )k\ +2k il Vikj -8ttV \{T ij + r ij + pu l u j ) 

(JjL 

+ (V^ku-Vikj + Vikl)k a (255) 
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(253) then gives using (255) to express the first term, (67) to express the second 
term, and taking into account the expression (251) of B J . 

dB j 

— = HB 3 + 2k\B l + ViW 3 (256) 

(JjL 

But by the Bianchi identities: 

1 
2 



Vii2 iJ ' = -V*i2 = g ij VjR = 



since R depends only on t. 

Finally (70) follows from (252), (256) and (250) 
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